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1 Introduction. 



It is well known that the structure of renormalized field theory is determined by the 
renormalization group (RG) equations. One considerable problem then is the calculation 
of renormalization group functions. Usually one can compute only a few orders in per- 
turbation theory. Hence, we do not have a complete knowledge of these functions. Other 
methods, such as the large N expansion, usually allow us to compute critical exponents 
at leading order only because calculations become too complicated at higher orders. 

To overcome these difficulties other methods have been developed. One such approach 
uses critical conformal invariance at the fixed point of renormalization group for calcu- 
lation of critical exponents, which are values of RG functions at fixed point. With this 
method it is possible to derive the critical exponents at order Using then the exact 

dependence of RG functions at each order in perturbation theory one can reconstruct 
them at some high-loop approximation. 

In this paper we derive five loop approximation for RG functions (up to few unknown 
coefficients) in the Gross Neveu and the non-linear cr models. 



2 Preliminaries. 

The [/(A^)-symmetric massive Gross Neveu (GN) model describes a system of D 
dimensional Dirac spinors, its non-renormalized action being 

S= U^x + mo)^. + ^7ofeV^a) V2] (1) 



a = 1, . . . ,N and summation over repeated indices is implied (in what follows we consider 
the theory in Euclidean space and functional distribution is denoted by exp 5* without 
the minus in the exponent). It is believed to be multiplicatively renormalizable in the 
framework of 2 -|- e expansion and its renormalized action in dimension D = 2 -|- 2^ is 

Sr = J d^a; [Zom^M'' + ^,{Z,d + Z^m)^^ + gM-^' Z^{i,Ja? /2\ , (2) 

where g is a, dimensionless renormalized coupling constant, m is a renormalized mass, 
M is a renormalization mass and ^0,1,2,3 are renormalization constants. The countert- 
erms of model (0) (without vacuum one) are generated by the standard multiplicative 
renormalization procedure 

^ Z^ip , mo = mZm , go = gZgM''^" , (3) 

Zi = Z'^ ■, ^2 = Z'^Zm , ^3 = Z'^Zg . (4) 

RG functions 'jp (the anomalous dimensions of quantities F) and f3 function are defined 
by the relations 



7ir = In Zp {F = ij,m,g) , (3 = V^g = {2eg - -fg)g , V^ = 



■ (5) 

90, mo 
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The (3 function of model (|^) has an UV fixed point [3{g^) = with (7* ~ e. Critical 
dimensions Ap = dp + '^*p (canonical + anomalous) of various quantities F (fields and 
parameters) in the framework of 2 + e expansion are defined by the relations 

AF = dp + 7; (F = V-, m) , A. = l/u = -(3'{g,) , (6) 

where = ■yp{g^), A^- = l/u is the critical dimension of "temperature" r = g — g^ and 
dp are the corresponding canonical dimensions: d^ = D/2 — 1, dm = ^■ 

For practical calculations it is more convenient to use the following model: 

S = Jd^'x [^,{d + mo)^a - a^l2go + ■ (7) 

The generalization of this model which is multiplicatively renormalizable in the framework 
of 2 + £ expansion is 



S = Jd^'x [^,{d + mo)^a - cry 2 + v^a^,^^ + gMJ,:fl2 + h^a\ (8) 

with two independent bare charges g^Q and (710 = Its renormalized form is 

Sr = J d^'x [Zom^M^' + i^^Z^d + Z^m)^;^ - ^3^V2 + 

+Z^ViM-'a^,^a + Z^g^M-^'i^.^af/'^ + ^e^MV + ha] . (9) 
If we integrate model (|^) over field a we obtain model (^ with parameters: 

m = m + hvM~^ , g = gi + g2 , gi = v'^ (10) 

and additional constant trln([9 + m]/[d + m]) + (/i^/2) Jd^x to the action. Compar- 
ing counterterms in the models (^ and (|^) one can obtain the following relations for 
renormalization constants 



Zi = Zi , = 1 + 2g,{Zo + n/8ne)_, Z2 = Z, = Z2Z3 

vZe = 1 - Z3 = 2gi{Zo + n/STre) , gfiZg^^Zl + 5(2^5 = gZ^ 



where gi = v"^ and n is defined by (|l^). These relations allow to express 7 constants 
Zi{gi,g2) of model (H) through 4 constants Zi{g = gi + (72) of model (0). It is more 
convenient to calculate the renormalization constants in the framework of massless (m = 
0,h = 0) one-coupling model like (|^) with interaction ~ ipipa and additional four-fermion 
counterterm which breaks the multiplicative renormalizability. 

To generate a systematic 1/n expansion one usually uses the model (|^) with a scalar 
field a. As have been shown in the massless version of (P) possesses critical conformal 
invariance at the Z)-dimensional fixed point of renormalization group and the exponent rj 
have been calculated up to order with the conformal bootstrap method. 

In what follows we shall consider 1/n expansions in arbitrary dimension of space D 
and the following notations will be used: 

NtTl = n, D = 2^, (12) 



2 



where tr 1 is a trace of the unit matrix in the space of D-dimensional spinors. In the 
framework of the 1/n expansion for the model (|^ the canonical dimensions of the fields 
if = ip,a and parameter r are 

d^ = d^ = fi- 1/2 , d„ = l, dr = 2{fi-l). (13) 

The notation d[F] = dp means canonical dimension of quantity F. Let us parametrize 
the critical dimensions Ap = dp + 'jp (canonical + anomalous) as follows: 

A^ = d^ + r]/2 , A^ = d^-r]-2A , Ar = l/u = 2X (Aq = /i - 1) , (14) 

where rj and u are the usual notations for the critical exponents, the quantity 2A is a 
critical dimension of the vertex ipipa in the model (^. For this model the quantity is 
connected with the critical dimension A[cr^] of a composite operator F = o"^ by means of 
the following relation 

Ar = l/iy = 2\ = 2i2- A[a^] . (15) 

For the 1/n expansion of any exponent z its coefficients in powers of 1/n are denoted 
by Zk'. z = zq + zi/n + Z2/'n? + . . .. The first coefficients ?7i, Ai, Ai for the special case 
of dimension D = 2^ = 3 have been calculated in [|I|]. The coefficient rji for arbitrary 
dimension D = 2fi have been calculated in p, Ai in ||^ and Ai, ri2 in [^]. In coefficients 
A2, A2 and Tjs have been calculated. 

The expression for 773 contains function /(/i) which cannot be expressed explicitly 
through the gamma-function T{z) and its derivatives. The function I{fi) is expressed 
through the self-energy massless diagram showed in Fig.l on which a line with index 
a means a massless propagator \x — According to its dimension the diagram is 

a simple line with index fi — 1 + A multiplied by the coefficient n(yU, A) which is the 
"value" of the diagram. The quantity /(/i) is defined from n(/i, A) by the relation J(/i) = 
d\nll{fj,, A)/(iA|A=o (for A = the diagram is calculated exactly: n(/i, 0) = Stc'^^ {ip' if^ ~ 
l)-V^'(l))r(/x-l)r(2-/i)/r(2/i-2)). in [|| the value of J(/i) for the dimension = 2/i = 3 
and the first terms of e expansions around D = 4 are given: /(3/2) = 3'?/^"(l/2)/27r^+2 ln2, 
/(2 — e) = 0{e) and in |^ the e expansion around D = 2 was obtained 

„l,,)._1^2C(3)£_,„),3,13C|)£!^,„ ,,,, 

An alternative way of derivation of e expansion for /(/i) is presented in the Appendix. 

Using 1/n expansion of critical exponents and some additional information about the 
explicit n dependence of RG functions one can derive high loop approximation for RG 
functions. This additional information usually includes the knowledge of exact expressions 
of the RG functions for some special values of n. 



3 Five loop approximation for the Gross Neveu and 
the a models. 

As have been noted in the case n = 1 {N = 1/2) corresponds to Majorana spinor 
and, hence, there are no four-fermion interaction term in (||), (|^), (H), (|^) and we have 
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for model (Q) 

Zo = 0, Zi = Z2 = 1 (17) 

and Z3 is undefined as coefficient of zero. 

Note that though the interaction {ijjipy does not exists at n = 1 but aipip does and 
for (|^) we have 

Zi = 1 , Z2 = = -^4 = 1 + — ! Zq = —— — Z2 , Zq = - — Z2 (18) 

Ane ?)2Txe Aire 

but again we can say anything about the value of at n = 1. 

From above considerations we conclude that for the model (^ ■j^(g)\n=i = 'Jm{g)\n=i = 
and, hence, these RG functions should contain the factor [n — 1) at each order of 
perturbation theory. But we cannot say the same about the (3 function. 

The case n = 2 corresponds to = 1, D = 2. As have been noted very early the 
U{1) Thirring and = 1 Gross Neveu model are equivalent in the two dimensions. The 
Thirring model in the two dimensions is known to be exactly soluble and it is easy to 
derive the RG functions (e = 0, A^ = 1): 

/?(«) = 0, 7^(«) = ^-^, u^j-. (19) 

But from it is known the following three loop approximation for RG functions in the 
model d): 

7^ = u\n - 1) - u\n - l){n - 2) + 0{u^) 

(3{u) = -2u\n - 2) + Au^n - 2) + 2u\n -2){n-7) + 0{u^) ^ ' 

setting n = 2 in it we obtain 

/?U=2 = 0,7v.U=2 = n' + O(M^) . (21) 

From the above results we can say that in the framework of 2 + e expansion the Thirring 
and the Gross Neveu models are not completely equivalent. Hence, there are no reason 
to believe that higher loop contributions of 7^ should contain the factor n — 2. 

Using now this information and results of it is possible to reconstruct RG functions 
of the Gross Neveu model in five loop approximation (up to some unknown coefficients) by 
expanding the critical exponents in both e and 1/n and comparing then the corresponding 
coefficients. The result is the following: 



7^(u) = u^{n-l) -u^{n-l){n-2) +u^{n-l)[n'^ -7n + 7] + 
W^^[3n=^(-2C(3) - 1) + 2n2(-3C(3) + 17) 

+12n(-llC(3)-10) + r75o] (22) 
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- 2)[-n^ - r2(66C(3) + 19) + Ao] + 

+ \u^{n - 2)[3n3(-2C(3) + 1) + n^(261C(4) - 168C(3) - 83) + 

+nP5i + /?5o] (23) 



7^(m) = 2u{n-l) ~2u'^{n-l) -2u^{n-2){2n-3) + 
in — \\ 

+2u^ ^ q n'(-6C(3) + 5) + n(90C(3) + 23) + m4o] + 
in — V] 

^ [3ri^(-6C(3) - 5) + n2(-81C(4) + 264C(3) + 169) + 



+nm5i + m5o] (24) 



Where 



7m = -7M=2e-7a (25) 

is the anomalous dimension of the mass and ^^so, /540 5 /^si, /^so, "^40, "^si, "^so are unknown 
coefficients. In 0] the expression (n — 2)(n— 5) = n^ — 7n+10 for the four loop contribution 
of 7^ have been received instead ri^ — In + 7. This result have been derived in using 
expressions for r/i, 772 and due to the supposition that multiplier {n~2) should be extracted. 
We checked our result (7 instead 10) by direct calculation of 30 relevant four loop diagrams 
in the MS scheme. So the supposition of extracting of multiplier n — 2 to all orders in 7^ 
RG function is wrong. 

Using a similar method and results of it is possible to derive the following expres- 
sions for RG functions up to five loops in the non-linear a model: 

'^2 2 8 

(n — 2) 

+M^^^^[-3n' + 2n + 14 + 3(3 -n){n + 5)C(3)] + 
(n — 2) 

+u'^——l[Qn'{3C{A) + 4C(3) + 2) - n'{5AC{3) + 5) 

-n(315C(4) + 1764C(3) + 307) + a] (26) 

P = 2eu-u\n-2)-u%n~2)-u'^^^^—^^^^-^- 
_^5 (^^-2) p^2 ^ 22n - 34 + 18(n - 3)C(3)] + 

^^6 (!1_^^3^3^2C(3) - 1) + r22(81C(4) + 144C(3) - 11) + bn + c] (27) 
96 

where a, b, c are unknown coefficients. 
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Appendix 

In this appendix we give the main steps of derivation of recurrent relation for function 
/(/i). This method allow to obtain the expansion (|l^) without the procedure of summing 
infinite series of the Gegenbauer polinomials The definition of /(/i) is 



A=0 



-2{m-1+A) 



^12 



where (see Fig.l) 

d X3ci X4 ^]^3 "^14 ^^34 "^23 ^24 ~~ "^(Z^; ' 

Using transformation from the function n(/i, A) can be transformed to 

n(/x,A) = c-F^(2/i-3 + A) , 



where 



and 



From 



J d X^d Xi Xi^ Xi^ X^^ X22; 2^24 — Ffj-io^) ' ^12 

r(^-i + A)r(^-i-A) 



(A.l) 

(A.2) 

(A.3) 
(A.4) 



1 + 0(A^) . 



r(/i-i)r(/i-i)r(i + A)r(i- A) 

it is obvious that the coefficient c is not contributed to /(/i) and we have 



(A.5) 



W= ^lnF^(2^-3 + A) 



(A.6) 



A=0 



Using the Gram determinant method |T^ one can obtain the following relation with 
shifted dimension of the space 



- F,ia) - F,{a + 1) + = ^^^^ + 2) 



(A.7) 



where Fj^^^a) is a sum of integrals which can be integrated using the chain rule. This 
sum is 



-An^[a^ + a(-3/i + 5) + /i^ - 5/i + 5] x 

^ r(/i - i)r(/i - i)r(/i - 2 - a)T{a + 3-2/1) 

~ r(a + 2)r(3/i - 3 - a) 



(A. 



On other hand, combining two recurrent relations [^] for the master diagram one can 
derive the following relation for F^(q;): 



F,{a + 1) = -jli^rizT^W + («) 



(A.9) 
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where is again a sum of easily integrated diagrams: 

2n^{3^i - 5 - 2«) r(/i - l)r(^ - l)r(/i - 1 - a)T{a + 4 - 2/i) 



^ ^ ' (2/i-3-a)(a + 2-/i)2 

From (|A.7|) and ( |A.9|) we get 

2(2^-l)(/i-a-2) 



r(a + l)r(3/i-4-a) 



F^+i(a + 2) + 



(A.IO) 
(A.ll) 



(/i — 2 — a) 



^ '^^ T{a + 2)T{3fi-3- a) ^ ' 



Substituting now a = 2/i — 3 + A and accounting that at A = (|A.2|) can be integrated 
exactly (see, for example, 0): 



r(2;U-2) 



(^V-i)-^'(i)) 



(A. 13) 



we derive the following relation for /(/i) 



B + 



6(/i-l) 



+ 



1 + 



+ + 



(A. 14) 



where 



= ^'(1) - ^'(^ - 1) , 5 = ^(2/i-2) + ^/'(2-/i)-^(/i-l)-^(l) . (A.15) 

Due to the property I{2 + e) = 0(e) f^, |^ one can obtain the expansion of /(I + e) 
{D = 2 + 2e) in e up to order e^: 



/(l + ^) = -| + ^e^-C(4)e^ + 0(.^). 



(A.16) 
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